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Abstract. A bosonized-field-theory representation of spin operators in the 
uniform-ficld-induced Tomonaga-Luttinger-liquid (TLL) phase in spin-1 Haldane 
chains is formulated by means of the non-Abelian (Tsvelik's Majorana fermion 
theory) and the Abelian bosonizations and Furusaki-Zhang technique [Phys. Rev. 
B 60, 1175 (1999)]. It contains massive magnon fields as well as massless boson 
fields. From it, asymptotic forms of spin correlation functions in the TLL phase 
are completely determined. Applying the formula, we further discuss effects of 
perturbations (bond alternation, single-ion anisotropy terms, staggered fields, etc) 
for the TLL state, and the string order parameter. Throughout the paper, we 
often consider in some detail how the symmetry operations of the Haldane chains 
are translated in the low-energy effective theory. 
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1. Introduction 

For (quasi) one-dimensional (ID) magnets, sophisticated theoretical tools have been 
developed in the last decades; for example, bosonizations El E|, conformal 
field theory (CFT) 0], Bethe ansatz |5| E3, and several numerical methods (exact 
diagonalization, density- matrix renormalization group (DMRG), quantum Monte 
Carlo method (QMC), etc). Thanks to them, we have obtained a deep understanding 
of ID spin systems. Of course, they have succeeded in explaining a number of 
experimental results of quasi ID magnets. In this paper, we try to further develop 
a low-energy field theory description for spin-1 antiferromagnetic (AF) chains in a 
uniform magnetic field, which are a fundamental model in ID spin systems. 

Before discussing spin-1 AF chains, for comparison, let us briefly review the well- 
established low-energy effective theory (so-called Abelian bosonization) for the spin- 
1/2 XXZ chain 0IH], which is one of the well-investigated and realistic models in spin 
systems. The Hamiltonian is defined by 

"1 
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; (StS- +1 +h.c.) + A z S*S* +1 



where Sj is spin-1/2 operator on site j, J > is the exchange interaction, A 2 is the 
anisotropy parameter (A z = 1 case is the SU(2) Heisenberg chain), and H > is a 
uniform field. In a wide parameter regime including the zero-field line (1 < A z < 1 
and H = 0) |oJ|8], the low-energy physics is governed by a one-component Tomonaga- 
Luttinger liquid (TLL) phase, where the spin correlation functions are of an algebraical 
decay type. For this critical phase, a beautiful field theory description has been 
established as follows. The low-energy effective Hamiltonian is a massless Gaussian 
model, 



Wf z = I dx 



K(dJ) 2 + hd x ^) 2 
K 



(2) 



where x — j x ao (ao is the lattice constant), 4>(x) is the scalar bosonic field, 9{x) 
is the dual of 0, v is equal to the massless spinon velocity of the chain and K 
is the TLL parameter (K = 1 and K = 1/2 respectively correspond to the XY (free 
fermion) point A 2 = and the SU(2) Heisenberg one A z = 1) 0. In this effective 
theory framework, the spin operator Sj is approximated [31 El IH1 d] as 
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where M is the magnetization per site (Sj}, and A n and A* L are nonuniversal 
real constants. The magnetization M, the velocity v, and the TLL parameter 
K can be determined as a function of the parameters (J, A Z ,H), via the Bethe 
ansatz integral equations for the XXZ chain. If these three parameters are given, 
one can correctly know the asymptotic behavior of any spin correlation functions 
through the formula (j3Jl . Since the accurate values of A n and A z n have recently been 
estimated ^2 EI], it is also possible to obtain the amplitudes of the spin correlation 
functions. This effective theory, especially the bosonized spin formula (0, is very 
useful in analyzing not only the XXZ chain itself, but also several perturbation effects 
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(bond alternations, next-nearest-neighboring interactions, anisotropy terms, etc) for 
the XXZ chain and coupled spin- 1/2 chain systems. 

Now, we turn to spin-1 AF chains and their effective theories. In this paper, we 
mainly consider the following simple spin-1 AF Heisenberg chain, 

Hs=i = J E $ ' - H J2 SI , (4) 

j 3 

where Sj is "spin-1" operator on site j. In contrast to the spin- 1/2 case, the spin- 
1 AF chain J3J without a uniform field (H — 0) possesses a finite excitation gap 
(call Haldane gap) on its disordered ground state, which is described by the famous 
valence-bond-solid picture j^]. The low- lying excitation consists of a massive spin-1 
magnon triplet around wave number p — tt/oq- The spin correlation functions exhibit 
an exponential decay. This drastic difference between the spin- 1/2 and spin-1 cases 
was first predicted by Haldane in 1983 [T3J: he developed the nonlinear sigma model 
(NLSM) approach 11411151 17] for spin systems in order to explain the above low-energy 
physics of the spin-1 chain. After Haldane's conjecture, based on the non-Abelian 
bosonization approach to ID spin systems ^JE1E|> Tsvelik proposed a Majorana 
(real) fermion theory |19l HI l2l I2()|. which also has an ability to demonstrate the low- 
energy properties of the spin-1 AF chains. Both theories are now a standard field 
theory method of studying ID spin-1 AF systems |23i like the Abelian bosonization 
for spin- 1/2 chains. 

When a uniform field H is applied in the spin-1 SU(2) AF chain the low- lying 
magnon bands are split due to the Zeeman coupling. Provided that the bottom of 
one band crosses the energy level of the disordered ground state, a magnon condensed 
state and a finite magnetization occur HU 1241 125] . It is well known that this 
quantum phase transition is of a commensurate- incommensurate (C-IC) type [261 15]. 
and the resultant phase can be regarded as a TLL. Therefore, it is expected that, just 
like the formula (0), one can obtain a bosonized-spin expression in this TLL phase. 
However, as far as we know, any clear bosonization formulas of spin-1 operators in 
this magnon-condensed phase have not been established yet. Namely, such a formula 
has not been argued enough so far |27| . 

The main purpose of this paper is to derive a satisfactory bosonization formula 
for spin-1 operators in the above TLL phase. To this end, we start with Tsvelik's 
Majorana fermion theory. This theory makes it possible to treat the Zeeman term in 
a nonperturbative way. Moreover, it can carefully deal with the uniform component of 
spins (the part around wave number p ~ 0) as well as the staggered ones, in contrast 
to the NLSM approach. (We explain these properties of the fermion theory in later 
sections). To make the new bosonization formula, the bosonized-spin expression for 
the uniform-field-induced TLL phase in a two-leg spin- 1/2 AF ladder [2Hj and its 
derivation, which was performed by Furusaki and Zhang |29| . is very instructive 
and helpful. We will rely on these results. We will often discuss the symmetry 
correspondences between the original spin-1 AF chain Q and the effective field 
theory. Like the formula Q), our resulting formula would be powerful in studying 
various (quasi) ID spin-1 AF systems with magnons condensed. Recently, magnon- 
condensed states in quasi ID gapped magnets, including spin-1 compounds such as 
NDMAP [mCIHE2E3ElES], NDMAZ |H], TMNIN |S7J and NTENP [SHI EH], have 
been intensively investigated by some experiment groups. Our new formula hence is 
expected to be effective for understanding results of such experiments. 

The organization of the rest of the paper is as follows. In Sec. |2 we review 
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Tsvelik's Majorana fermion theory for the spin-1 AF chain without external fields. 
Section |3| is devoted to discussing the effective field theory for the uniform-field- 
induced TLL in spin-1 AF chains, which is based on the Majorana fermion theory. 
These two sections are needed as the preparation to construct a field-theory formula 
for spin operators in the TLL phase, although several parts of them have been already 
discussed in literature [Tl l2l 1191 1201 140) . In Sec. 01 which is the main content of this 
paper, employing the contents of Secs.|21and|ni we derive the formula of spin operators. 
Moreover, using it, we estimate the asymptotic forms of spin correlation functions in 
the TLL phase. Section[S]provides easy applications of our field-theory formula of spin 
operators. We consider the nonlocal string order parameter |411 140j . effects of some 
perturbation terms (bond alternation, a few anisotropy terms, etc) for the TLL state, 
and the magnon-decay process caused by anisotropic perturbations. We summarize 
the results of this paper in Sec. El Appendices may prove useful when reading the 
main text. 



2. Spin-1 AF chain without external fields 

In this section, we review the non-Abelian bosonization (Tsvelik's Majorana fermion 
theory) method for the spin-1 Heisenberg AF chain with H = 0, especially focusing 
on the symmetries of the chain. 



2.1. WZNW model and Majorana fermion theory 
First, we define the spin-1 bilinear-biquadratic chain, 

fip = J^2 [Sj ■ S j+1 + P{S 3 ■ S j+l f\ . (5) 

3 

Several studies 021 reveal that the system belongs to the "Haldane phase" in the 
regime \[3\ < 1, and the points (3 = ±1 correspond to a quantum criticality. In the 
Haldane phase, the ground state is disordered (namely, all symmetries are conserved) 
and the lowest excitations have a finite Haldane gap on it. It is shown [161 1431 134] that 
the low-energy physics of the critical point (3 = —1 is captured by an SU(2) level-2 
Wess-Zumino-Novikov-Witten (WZNW) model [1I2J, a model of CFT with central 
charge c = 3/2, plus irrelevant perturbations preserving the spin SU(2) symmetry. 
This WZNW model possesses a Majorana (real) fermion description, i.e., the WZNW 
model with c = 3/2 is equivalent to three copies of massless real fermions, each 
of which is regarded as a critical Ising system (a minimal model of CFT) with 
c = 1/2 ^3|2HinEl- Using the fermion (or Ising) picture, one can represent the 
Hamiltonian of the WZNW model as follows: 



« ~ 1 O Q 



9=1,2,3 



where Cl(r)( x ) = ^l(r) ( x ) ^ s tne ^ (right) moving field of the qth real fermion with 
scaling dimension A s = 1/2, and vo is the "light" velocity of the fermions. We here 
normalize the fermions through the anticommutation relation {£$(x), (y)} = 

d u yS qjq i5(x — y), where the delta function stands for 5 x ^ y /a [a is a constant of 0(oq)]. 
At the WZNW point j3 = — 1, the spin operator is approximated as the following sum 
of the uniform and the staggered components |lfil 1191 ^ 12] > 

S]/a Q ~J q (x) + (-iyC Ni(x), (q =1,2,3 or x,y,z), (7) 
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where 

Ji{x) = Jl(x) + J%(x), Jl {R) = -i?%k)%&v ^ + 3 = «)' 
N 1 = f^a 2 ^ 3 

N 2 = a 1 fi 2 f x 3 , (8) 
N 3 = /iV 2 ^ 3 

and Co is a nonuniversal constant. The quantity Jl(r) IS equivalent to the left (right) 
component of the SU(2) current in the WZNW model with A s = 1, and a q (/j, q ) is 
the gth-Ising order (disorder) field with A s — 2/16. 

According to Tsvelik [T^l UHl HI 12], transferring from the WZNW point /3 = -1 
to the Heisenberg one {3 = corresponds to adding two SU(2)-symmetric (see the 
next subsection) perturbation terms, a fermion mass term and a marginally irrelevant 
current-current interaction, to the WZNW model ©. Namely, he proposed that the 
low-energy and long-distance properties of the spin-1 Heisenberg chain J3J are govern 
by the effective Hamiltonian, 



(cm - &d x &) + mie L e R \j q L j q R 



, (9) 



where we may set m > and A > 0, which guarantees the irrelevancy of the 
A term. From the viewpoint of the spin-1 AF chain, fermions R stand for 
the low-lying magnon triplet around wave number p = tt/oq, the mass parameter 
m corresponds to the Haldane gap, and the A term represents the inter-magnon 
interaction. Renormalizing the A interaction effect into the kinetic and the mass 
terms, or simply neglecting the A term, one can determine the values of vq and m from 
the numerically-estimated magnon dispersion of the spin-1 Heisenberg chain; such a 
procedure provides vo s» 2.49Jao and m « 0.41 J ^HlEEj- One should interpret that 
the inequality m > corresponds to the disordered phase in the Ising-model picture, 
in which (/z 9 ) ^ and (a q ) — 0. For the disordered phase, Ising-held correlation 
functions |17| |23 EH are shown to be 

(<7«(xK(0)) « -j== (1 + Odx/ar 1 )) + 0(e~ 3 ^), 



+ 0(e- 4 W/S°), (10) 

at a long distance, x 3> £ c = va/m |48| . Here, B\ is a nonuniversal constant, and 
£ c s» 6oq is the correlation length. If we assume that the field-theory formula (J7J) 
is valid even at the Heisenberg point /3 = 0, that formula and the result IjlOjl lead 
to the following asymptotic behavior of the spin correlation functions of the spin-1 
Heisenberg chain: 

(S q S q ) « f IV B2 _j_ ^3 -2|x| /g, , . fill 

( J o) ( } /M7£ N/6. ' ( } 

where we used Eqs. (|B.3|) and (|B.4(I to derive the second term, and B2.3 are 
nonuniversal constants. The exponential parts in Eq. Ijlljl are consistent with those 
predicted by the NLSM approach , but the prefactor (cx 1 /x) of the second uniform 
term differs from that of the NLSM, which predicts a prefactor proportional to 1/x 2 . 
From the calculation of DMRG (Hj, we may conclude Bi ~ 0.25 and B3 <C Bi- 
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2.2. Symmetries of the spin-1 AF Heisenberg chain 

The effective Hamiltonian © and the field-theory expression of the spin 10 tell us 
how the symmetries of the spin-1 Heisenberg chain are translated in the field theory 
language |161 1181 |2"U| . Here, we consider four kinds of symmetry: the spin rotation 
Sj — > T^-Sj (1Z is an SO(3) matrix), the one-site translation Sj — > Sj+i, the time 
reversal Sj — ► — Sj, and the site-parity transformation Sj — > S_j. Note that the 
link-parity symmetry is the combination of the one-site translational and site-parity 
symmetries. 

(i) The spin rotation may be realized by 

(l,r -» K£ LtR , N - UN, (12) 

where C = t (il,£,l,Q) and N = t {N 1 ,N 2 , N 3 ). The "rotation" of fermions induces 
J v -> TIJ V and J — > 72.J, where X = '(J*, J 2 , J 3 ) and J = '(J 1 , J 2 , J 3 ). The latter 
transformation in Eq. I|12|) could be regarded as 

a^TZS, (13) 

where <? = *(<r 2 , cr 1 , c 3 ). If both the mass and the interaction terms are absent (i.e., 
the system is completely critical), the chiral SO(3)xSO(3) symmetry (independent 
rotations of the left and right fields) appears. 

(ii) The one-site translation may correspond to 

Cl,r( x ) -> -bM* + «o), N(x) -> -N(x + oo). (14) 

Since the fermions R describe the magnon excitations around wave number p = ir/ao 
in the original chain (J3J, we should insert —1 in the first transformation of Eq. |Q - 
The second transformation can be realized by 

a(x) — > —a(x + ao) j 3{x) — > —a(x + ao) 



(15) 

p(x) — * fi(x + ao) ' " i \ fi(x) — > — p(x + ao) ' 

where fi = '(/x 1 , fx 2 , /i 3 ). Due to this transformation, the factor (— l)- 7 in front of the 
staggered component of the spin Q changes to (— 1) 3+1 , 

(iii) The time reversal may correspond to 

N^~N, *--». (16) 

The second mapping AT — > —AT could be interpreted as 

^ r , or <^ ^ ^ . 17) 

(iv) As a proper mapping of the continuous fields towards the site-parity 
transformation, we can propose 

I ff^+ftl < (18) 
{ (r{x) -> ±&,(-a0 

All the transformations (i)-(iv) leaves the effective Hamiltonian invariant. In 
other words, the symmetries (i)-(iv) allow the existence of the mass term and the 
current-current interaction one. 
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2.3. Dirac fermion and Abelian bosonization 

Besides the relationship between the WZNW model and three copies of massless 
real fermions, it is well known that two copies of the real fermions are equivalent 
to a massless Gaussian theory with TLL parameter 1 or a massless Dirac (complex) 
fermion, except for a few subtle aspects 0]. Namely, one can bosonize two copies of 
the real fermions. In this subsection, we apply such a bosonization procedure for the 
effective theory ©. The results are useful in the later sections. 
Let us define the following Dirac fermion, 

1>L = fye L +i&), 1>r = ^=(& + *&)- (19) 

which obey {i> v (x), Vv (y)} = and {ip v (x), V'Mj/)} = S u y6(x - y). Using this, we 
can rewrite the bilinear (free) part of the Hamiltonian J3J as follows: 

Hfroo = J dxiv (iJ^dxijjL - ^ndxtpjij + mi{4 L i , R - ip R ipi) 

= %>]+% 3 ]. (20) 

Applying the standard Abelian bosonization method (see |Appendix A| to the Dirac- 
fermion part Tify], we obtain the following sine-Gordon Hamiltonian, 

Hty] _> 7i[$] = f dx ^ [(d x Q) 2 + (d x <S>) 2 ] + — cos(V4^$), (21) 
J I not 

where $ is the scalar field, 9 is the dual of <£>, a ~ ao is the nonuniversal short- 
distance cut off, the relevant cos term bears the Haldane gap, and we simply dropped 
the Klein factors. The fields and (cr 1 ' 2 , A* 1 ' 2 ) in the spin operator (JJJ can also 

be rewritten in the Dirac-fermion or the boson languages. The uniform component in 
Eq. is 

Jl + Jr. - -^(4 ~ 4>l)& + 4>R)& 

= -}= [( L e L sin(v^(<i> + 9)) - ( R e R sm(M® - 9))] , 



ji + j 2 r = ^=(vi + «i + -^A + iPrKr. 



[CUl cos(V^($ + 9)) + Ck& cos(V^(<i> - 9))] , 

4 + Jr. = - ■ 4^ ■ - ■ 4 R ^R ■= -A=d*$, (22) 

It; 



where is the Klein factor for the fermion ip v , and we assumed {Cv,^t'} = 0. The 
symbol : : denotes a normal-ordered product; : ij^lipv '■= 4>t>'4>v — (V'JVv)- From the 
formula (|A.6|I . the staggered component of the spin [S3j is bosonized as 

N 1 = n 3 sm(y^<d), 
N 2 = ^ 3 cos(^F9), 

N 3 = a 3 cos(V^^). (23) 
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Employing Eqs. (|19|l - i|23|) an d Abelian bosonization formulas (see |Appendix A} , 

let us consider how symmetry operations of the original spin-1 chain Q are represented 
within the above Dirac-fermion and boson framework, as in the last subsection, 
(i) Let us define the following U(l) rotation around the spin z axis, 

(24) 





where 7 is a real number. One can easily find that this corresponds to 

i>L,R -» e l > L ,fl, ®{x) -» Q(x) ~ 1=. (25) 

v 7r 

This transformation is nothing but the U(l) symmetry that is explicitly preserved via 
the Abelian bosonization procedure. 

(ii) The one-site translation Sj — > Sj+i could be realized as 

iPl,r(x) — ► —iPl,r(x + ao) ( 9(x) — > Q(x + ao) + y/n 

£UM ^ -£U0£ + ao) ' 1 $(x) -> $(x + a ) 

^r"" 3 /^ ^ • (26) 

The second and third transformations change the staggered factor (— l) 3 in Eq. (J7J 
into (— 1) J+1 . These two correspond to the first proposal in Eq. (|f 5fl . 

(iii) The time reversal & — ► — could be interpreted as 



R,L 

3 



1pL,R *-> i>l 

a 3 -> -a 3 
/i 3 —> -^ 1 




(27) 



The second and third mappings are equivalent to the second proposal in Eq. Q17[l. 
(iv) The site-parity transformation Sj — * S-j might correspond to [HJ 

■0ij(a:) -> ±ifi L (-x) 



aiW^TtM-x) 

a 3 (x) — » cr 3 (— x) 
At 3 (x) — > ^i 3 (— x) 




(28) 



3. Effective field theory for uniform-field-induced TLL in spin-1 AF chain 

Based on the Majorana fermion theory explained in the last section, we consider the 
uniform- field effects in the spin-1 AF chain 0}, especially the uniform-field- induced 
critical phase. 

From the formula J7}, the Zccman term in the chain 0} is approximated as the 
following fermion bilinear form, 

- H E s i ~ / dx w&zi + e R e R )- (29) 
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Therefore, it is possible to nonperturbatively treat the Zeeman term, together with 
the free part Q2U[) . within the field theory scheme. This is a significant advantage of 
the Majorana fermion theory. Take notice here that the Zeeman coupling reduces the 
spin SU(2) symmetry to the U(l) one around the spin z axis, and destroys the time- 
reversal symmetry; the remaining symmetries are the U(l), the one-site translation, 
and the site-parity ones. The right-hand side in Eq. 129(1 is indeed not invariant under 
the time reversal and any spin rotations except for the U(l) one. 

Supposing the inter-magnon interaction may be negligible at the starting point, 
the effective Hamiltonian for the spin-1 AF chain Q with a finite H is 



Tif= I dx 



iv (^ L d x t/jL - tpjtdxipR 



+ 5«o (ddsfi - e|0x£|) + mi&&] . (30) 

In order to diagonalize this Hamiltonian, and see its band structure, we introduce the 
Fourier transformations of fermion fields as follows: 

= ^Y,e lkx iUk), H(k) = -±=J dxe-**l*{x), 

i> v {x) = -±= J2 e lkx Mk), Mk) = -^=Jdx e~ tkx Mx), (31) 

where L — Nag is the system size (N : total number of site), and k (\k\ < A : A ~ a^ 1 
is the ultraviolet cut off) is the wave number: a related quantity p = k + n/ao is 
interpreted as the wave number of the original lattice system Since Q is real, 
Q (k) = k) holds. New fermions t$ and ip u obey the anticommutation relations 
= Su,wSk,-v, bi> v {k),i>l,{k')} = 6 v y5 k ,«, and {^(fc),^(fc')} = 0. 
Substituting Eq. (|31|l into the Hamiltonian (|30|l . we obtain 

nf= £ H^ k + £ %Mis k) (32) 

|fc|<A 0<fe<A 

where * fc = *(^,^l), E k = s (ffj,f£), and 

Mt = ( kV0 + H ~™ A M{=( S -£ ) . (33) 



fe Y im —kvo + H J ' k \ ' lra ~kV\ 
Subsequently, performing the following Bogoliubov transformation 



where 



HS V ), (34) 



V+(k) J ' ' V U-k) 



(eo(fc) - kvo)- 1 / 2 (co(fc) + fc« r 1/2 



v /2i^(fc) V *( £ o(fc) ~ fct>o) 1/2 /™ -*(co(fc) + fcv ) 1/2 /™ 



eo(fc) = ^Wvl+m\ (35) 
one can finally obtain the diagonalized Hamiltonian, 

Kh= E [eAkW-(k)v-(k)+eo(k)fjUk)fj (k) 
\k\<A 



(k)fil(k)fi + (k) , (36) 
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where e T (k) = e (/c) ±H. While three bands e±o(k) are degenerate at H = 0, a finite 
field H mixes the first and third fermions and £ 3 , and leaves these bands split. The 
dispersions e± o(fe) indicate that the Majorana fermion theory correctly reproduces the 
Zeeman splitting. Fields rf^ and fj_ may respectively be regarded as S z = 1,0 and 
= — 1 magnon creation operators. Here, let us consider how the symmetry operations 
are represented in the picture of the magnon fields ?7±,o- Using the transformations 
<l25l) and l (25jl and the relationship between (^„,£j)) and f}±fl [Eqs. i|3*T|) and iJSU], 
one finds that the U(l) rotation along the spin z axis and the one-site translation 
respectively correspond to 

fj ± ^e^fj ± , (37) 
*°a*.rn J -e 2feao ^o(fc) (fe>0) 



If we neglect the fc-dependence of or consider only fermion fields around k = 0, we 
could interpret that the site-parity transformation in Eq. I|28[) is realized by 

The Hamiltonian l|36|l is clearly invariant under these transformations. One has to 
keep in mind that the symmetry operation i|39J) is valid only around k = 0. 

As 7? exceeds the Haldane gap m, the band e+(fc) crosses the zero energy line. 
This just corresponds to the C-IC transition of the spin-1 AF chain (gj [221 1%3"1 1211 125] . 
Then the system enters in a critical TLL phase with the magnon fy + condensed. The 
two remaining bands e_ i o(fc) are still massive and completely empty of magnons: 
(/i 3 ) 7^ and (cr 3 ) = still hold. Focusing on the low-energy physics in this case of 
H > m, we can approximate the S z — 1 magnon field fj + as a new Dirac fermion with 
a linear dispersion as in Fig.Q] The left and right movers L(x) and R(x) of the Dirac 
fermion are defined as 

L(x) = -L J2 e lkx m, R(x) = 4f E e <te fl(fc), 

Viy |fc|<A' Viy |fc|<A' 

L{k) = fi+(k + k F ), R{k)=f l +{k-k F ), (40) 

where the Fermi wave number kp = \/ H 2 — m 2 /vq is determined from e + (kp) = 0, 
and the cut off A' should be much smaller than A. Using the Dirac fermion, one can 
rewrite the effective Hamiltonian (13611 in the case of H > m, in the real space, as 
follows: 



^ff>m — I dx 



IVp 



(tfd x L - R ] d x R) 



■ri 



-^d 2 x +m + H + 0(d 4 x ))v- 



2 



= H[L, R}+H[e\ +H[r,-), (41) 
where vp = de+/dk\k=k F = \/H 2 — m 2 vo/H is the Fermi velocity, and r)-(x) = 
-j= e tkx Tj- Take notice that the linearized dispersion and the definition of L 
and R become less reliable as the field H becomes closer to the lower critical value m 
{vp -> 0). 
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Figure 1. Spin-1 magnon band structure at H > m. 



So far we have omitted the inter-magnon interaction terms in this section. They 
would yield interactions among L, R, Q and ry_ (see the next section). A better 
method of dealing with such terms is the Abelian bosonization. Through it, the 
Dirac-fcrmion part in the effective theory 1411) is mapped to the following Gaussian 
model, 



H[L,R] -*H[4>] = / cfo| 



K{d x 6) 2 + hd x 



(42) 



where <f> is a scalar field, 9 is the dual of <j), v is the renormalized Fermi velocity, 
and K is the TLL parameter: although K = 1 and v — vp at the stage of Eq. I|41|) . 
both parameters are subject to the renormalization due to irrelevant interaction effects 
neglected here. The Hamiltonian (|42[) just describes the low-energy physics of the TLL 
phase in the chain 

4. Field-theory representation of spin operators in H > m 

Making use of the contents in the preceding sections, we construct a field-theory 
representation of spin operators in the uniform-field- induced TLL phase in the spin-1 
AF chain (gj. 

4--1- Spin uniform component and Magnetization 

First, we consider the uniform component of the spin, J = Jj, + Jr in Eq. (JJJ. In the 
low-energy limit, Eqs. I|34|l and l|40(l allow us to approximate the "old" Dirac fermion 
as follows: 

xp R (x) » -=ri-(x) + e lkFX U+L(x) + e - lkFX U^R(x), 
V2 

i/) L (x) » -^v-(x) - ie %kFX U-L(x) - ie~ ikFX U+R{x), (43) 
where the prefactor l[i]/V2 of the first term is the (1,1) [(2,1)] component of the 



matrix Uk=o, and U± = m/ y 2H(H ± V H 2 — m? ) is the (1, 1) component of U T k F - 

Note that U± have properties (i) U± l/y/2 at H -> m + 0, (ii) U\ + Ul = 1, 
and (iii) U+U- — jjj. From this relationship between and (L,R,rj-) [see also 
Eqs. I|37|) - (|39|) ]. we know how the latter three fields are transformed by symmetry 
operations. The U(l) rotation in Eq. I)25H corresponds to 
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The one-site translation in Eq. I|26|) is mapped to 

L(x) -> -e ikFao L(x + a ) 



rj-(x) -> -r)-(x + a ), 



R(x) -> -e- lkFa "R(x + ao) 



(45) 



Tipn(-x) and ip R (x) 



The site-parity transformation in Eq. I|28|l . 

±iPl(—x), can be reproduced by 

r,_(x)-.±^_(-x), | i?(a;) _ >TjL( _ x) • (46) 

One can confirm that the Hamiltonian l|41l) is invariant under these transformations. 

Using the relation l|43(l . we can rewrite fields £*' 3 in J^/j in terms of L, R, and 
r/_. As a result, the currents Jl.r are reexpressed as 



~{V- + V-) + -^U-{er ikpx L^ + e lkFX L) 
2 v2 



+ -^-U+(e ikFX R^ + e- lkFX R) 
l(r ? l- 7? _)--L C /_(e-^Lt 
1 



£ 3 



(47a) 



-^U+(e ikFX R^ -e- lkFX R) 



£ 3 



(476) 



V2 



rf_L + h.c 



V2 



ikpx f 



r/li? + h.c 



- C/+Z7_ [e- 2ifel?x i t ii + h.c] , 

- T)-) + -^=U + (e- ik ™tf - e lkFX L) 



1 



+ —U-(e tkFX R^ - e- ikFX R) 
v2 



£ 3 



(47c) 



{Aid) 



+ n-) + ^U+{e- lkFX L^ + e lkFX L) 
2 v2 



ikpx j^f 



R) 



3 



(47e) 



V2 



e ikFX nlL + h.c 



V2 



— ikpx t 



7/1-/? + h.c 



[/+[/_ [e-^'i+JZ + h.c] , 



It is easily found that the above Jl.r are appropriately transformed for the symmetry 
operations f4*4*|) . (|4l))) . (J2SJ|, (|4^|l and (J2BJ. In addition, the Hermitian nature of 

all current operators is preserved in Eq. (47). Even if, instead of Eq. (|43l) . we use 
precise representations of continuous fields in the Fourier space such as Eqs. (1311) and 
<|4(Jll . we can finally arrive at the same expression as Eq. (47). 

The Fermi wave number fc^ is related with the uniform magnetization M — (Sj). 
It is given by the formula (47) or the Fourier-space representation of J 3 = j£ + J R . 
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mX> s *) = Jff dx<J * ix)) 









fdx\ 


V NJ 





\k\<k F 



(Lit - tfL) + {RR ] - R)R) + (v-V-) 

kpCLQ 



(48) 



From the derivation process of Eq. Q48fl. we also find that the first three terms in 
Eqs. I|47e|l and Q47j) should respectively be regarded as 



zap 2 



2 



- ( 77-77! - rjlr)- 



^± (LLt - LtL) + ^ (JIRt _ fltij) 



-Ul : L f L : -t/ 2 : B*R : +— - V 
+ " 2a 2 ' 



(49) 



Utilizing the results (47) - (|49|l . and applying the Abelian bosonization to the 
Dirac fermion (L, R), we can straightforwardly lead to the following partially bosonized 
currents Jl,r- 



Jlix) 



^(V- + + l ^— ^L cos(V47t0l ~nMx/a ) 

2 \ 7r<y 



H — + Kg cos(v / 47r(ftfl — irMx/ao) 
Vita' 



f 3 

SL> 



TiiV- - V-) 1=f Kl sin(V47r0L - irMx/ao) 

2 yira 



iU+ 



'ira 



kji sin(v47T(/)fl — nMx/ao) 



£ 3 

Si! 



f/ 2 



y7r zap 2 



/7T 

?7_ 



V47T0 7 

u+ 



7?t ^^g-^v^iAi+'iTrAfx/ao _|_ 

^^iv^A-mMs/ao + h.c 



V47I 

— i — — klkr sin(V^7r0 — 2irMx/ao), 

na 

Jr( x ) ~ 7t( ? 7- ^ ? 7-) + r— KL sin(V47T(/)L - irMx/ao) 

L2 V7TQ! 



iZ7_ 



sin(v 4^0^ — nMx/ao) 



(rjl_ + 7?_) + —j=^=.kl cos(\/47r0 L - nMx/ao) 
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iU_ 



kr cos(v Att<Ph — irMx/ao) 



JUx) « - 



u 2 


L ~ 


U 2 


u+ 


rj_ 


K L e~ 


\f Attol' 


U- 


V- 


K R e %s 


V Ana' 


.U+U- 







d x 4>R + 



M 
2a 



1 



ira 



(50) 



where a' ~ 1/A' is the short-distance cut off, and Kl(r) IS the Klein factor for the 
field L(R). Therefore, the bosonized uniform components of the spin are written as 



J 3 (x) rs — \=d x (j> — r)\_r)- — iCqRl^r, sin (V47T0 — 2-KMx/ao) 

a V 71 " 

+ Cf ^lK L e- i ^ F(<!i+6 ' )+4 ' rM:E/ao + h.c) 
- C? (r 1 l KR e- l ^- e) -™ Mx/ao + h.c) + • • • , 
J+{x) = J 1 (x)+iJ 2 (x) 



(51a) 



it]- + e 



+ iC 2 R R e i ^- vMx/aa + ■ ■ 



f 3 



+ C 1 K R e l ^ 4 '-' rMx/ao + 



(516) 



where C n and C^ are nonuniversal constants, which contain the effects of irrelevant 
terms neglected at the stage (|41f) . When the irrelevant terms are not taken into 
account, constants C§, Cf, Ci, and C 2 are respectively proportional to U+U-, 
U- — U+, U-, and U+. Therefore, it is inferred that Cf — > at H — ► to + 0. In 
Eqs. I|5()[l and (51), we used the formula (| A. 3|) . If, instead of it, we use another 
formula (| A.4|> . we can add more irrelevant terms to Eqs. (|50(l and (51). Namely, we 
may perform the replacement, 



a ±iy/T?($+B)^iirMx/a 



^±i^/7r((fi—6)^iTrMx/ao 



n=0 



where C n is a nonuniversal constant. Equation (51) is a main result in this paper. 
We expect that replacing the Klein factors in Eq. (51) with a constant is admitted in 
most situations. (A typical situation is when one investigates asymptotic features of 
spin correlation functions. See Sec. l4.4h 

Finally, using the results (47V 15U|) . let us investigate how the current-current 
interaction A term modifies the parameter K of the effective theory (14211 . Since the 
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SU(2) spin-rotation symmetry is reduced to the U(l) one in the case with a finite 
H, the term —XJlJr might be modified to —X±(J^J^ + J\J^} — Ay J\J\, in which 
we assume that both Aj^ and Ay are positive. Equations (47) and l|49|l lead to the 
following representation of the current-current interaction |52j : 

Jl(x+)J%(x) « {U 2 _ : R)R : +U% : L^L :) (U\ : R ] R : +Ul : L ] L :) 

+ U\ Ul (L^R(x+)R^L{x) + R^L{x+)L^R{x)) 

+ XP_ 

r\'_r] 

(L(x+)tf(x) + R(x+)ltf(x)) + h.c. 



+ - nln- (■ L^L: + : R^ R 

U+U- 
2 

U 2 _ -Ul 



(rjl_ {x + )i]l_ (x)LR + h.c.J 



2 

+ oscillating terms + others, (53a) 

Ji{x+)J^x) + Jl{x + )J 2 R {x) « 

iC/+C/_ (i f L - Lit + - flflt) ^£^3 

+ oscillating terms -I- others, (53 6) 

where we used x+ = x + S (5 is a small parameter of O(ao)), "oscillating terms" have a 
factor e ±m7rMx / a o ( n: integer), which are irrelevant except for the case that M is close 
to a special commensurate value, and "others" are constructed by only massive fields 
(£„,77_). The Hamiltonian l|41l) plus the above current-current interaction is regarded 
as a more accurate low-energy effective theory for the uniform-field-driven TLL phase. 
Integrating out the part of massive fields , ry_ ) in the effective theory [S3] , we obtain 
a Hamiltonian for the interacting Dirac fermion (L, R) that just corresponds to the 
Gaussian theory (|42|l . In the Gaussian theory framework, the first and second terms 
of Eq. (|53a|) directly contribute to the renormalization of the TLL parameter K and 
the velocity v as follows. The bosonized forms of these two terms are written as 

(U 2 _ : R^R : +U% : L^L :) (U% : i? f i? : +U 2 _ : L^L :) 

1 ~ l ^[d x <t>d x 6-d x 6d x <j>], (54a) 



4tt 

{/it/ 2 . (L ] R{x + )R ] L{x) + R^L(x+)iSR(x)) 



U\U 2 _ 
{2-Ka'f 
m 



(d^) 2 + const + • • • , (546) 



where we used U\ + Ul = 1 and U+U- = jjj, and assumed 5 — a' . The final term 
in Eq. H54af> could be negligible or be replaced with a constant: such a procedure 
is supported by the equal-time commutation relation [<j)(x),6(y)] = —|sgn(x — y). 
Consequently, two terms in Eq. (54) lead to an additional boson kinetic term, 

- ~ ^Ax)(^) 2 + ^(Ul Ul)\d x 9f, (55) 
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for the K = 1 Gaussian theory 1|42[1 [53| . This term is inclined to make the parameter 
K increase, i.e., K > 1, provided that An is the same order as Aj_. The relation K > 1 is 
consistent with the predictions of the numerical calculation in Ref. [25] and the NLSM 
approach ,24 . (It is shown in Ref. (25] that K is between 1.0 and 1.5 in the field- 
induced TLL state.) Besides Eq. (54), other terms in Eq. (53) also yields a correction 
of K and v, but it is expected that the property K > 1 is maintained. From Eq. (53), 
one also finds that the current-current interaction does not generate any relevant or 
marginal vertex operators such as cos(V4tt0). It means that the massless TLL phase 
survives even when the A interaction is taken into account. The third normal-ordered 
term in Eq. I|53all provides only a small correction of the magnetization M, and can 
be absorbed into the Gaussian theory via <f)(x) — » 4>{x) — const x x. 

4-2. Spin staggered component 

Compared with the formula (51) for the uniform component of the spin, the way 
of leading to a field-theory formula for the staggered component N(x) is not 
straightforward, because we have to express the vertex operators e ±l ^ e and e ±4v/ ^* 
in Eq. i|23|) in terms of another boson ((/>, 9) language and the magnon 7/_ one. To this 
end, we apply the discussion in Ref. |29|. 

First, using some results in the previous sections, let us naively speculate the 
relation between the "old" boson ($,0) and the "new" one (<j),0). The expression of 
J 3 {x) in Eqs. I|22l) and (|51a|l indicates the following relation between $ and 0, 

- -?j=d x ®(x) w - -=d x 4>(x) + M/a 

V7T \/7T 



- iCoK L n R sin(V47n/> - 2TrMx/a Q ) H . (56) 

Performing the integration with respect to x in Eq. I|56|) , we obtain 

» 4>(x) - y/nMx/ao + ■ ■ ■ , (57) 

where we dropped the sin term because it oscillates and then vanishes in the 
integration. From this, we can expect that 

e ±iv^F$ ^ e ±i(- v /7r0-7rMa:/oo) ) (53) 

holds within a rough approximation. We next consider dual fields and 0. The U(l) 
spin rotation in Eq. 14411 could be realized by 

0(x) -> 9{x) - -1=, (59) 

in the boson language. The comparison between this and Eq. i|25|) implies the relation, 
e ±*V¥e _ e ±V5» (60) 

In order to raise the validity of speculations l|58() and l|60[) . and obtain a more 
appropriate relationship between e ±*\/?f*(©) anc j ((f> 7 Q 7 r]-), a standard bosonization 
formula, 

</>L = -^=e-^(*+e) 5 ^ = Cfl^HM*-^ (61) 
\J2-kol \/2na 

is available. From this formula, we rewrite vertex operators e~ iv/ ^* and e~ l ^ & as 

e" 1 ^* = V2^ [ve^ e C L ^L + (1 - 2?)e- 4 ^ e CiiV4] > (62a) 
e-^ e = V2^ \ve l ^<; L ^L + (1 - V)e-^( R <p R ] , (626) 
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where T> is an arbitrary constant |55) . Then, employing approximated results l|58|) . 
<|6U[I . (I43|> . and the bosonization formula L,R ~ e Tiv^<t>L,B m the right-hand sides in 
Eq. (62), we can arrive at a desirable representation of e ±lv/ ^*( e ). This method was 
proposed in Ref. |2*5| . 

Following the above idea, we obtain 



V- 



U- {-KlKl + Cfl«fl) e~ l ^ + * Mx ' a ° 

+ C/+ (-*C £ «B + Cfl«L) e l ^-A/,/ao 

+ ^(^-e i ^ + C«r ? ie-^)+-- 1 



(63 a) 




(636) 



where we set V = 1/2, and is the Klein factor of the field L(R) |5B|. It is 

found that the naive expectations (|58|) and l|6U|l are qualitatively consistent with the 
result (63). From this formula, our target, the staggered component of the spin may 
be represented as 

N 3 (x) = cr 3 cos(\/7r$) 



« fi 3 7?i cos - 7rMa:/ao) 

+ Df / {(<a + C«)»?-e i ^ f9 +h.c.} 

A^+(a;) = TV 1 ^) + i7V 2 (a;) = i/i 3 e"^ e 

« ^e~^ e Di + D 2 cos (Vl^d - 2TvMx/a 



(64a) 



+ 



(646) 



where D„ ' are nonuniversal constants, which include effects of irrelevant terms, 
and we naively omitted some of Klein factors and k v . We used the formula 
L,R~ e Ti^<t>L, R m gqs. (63) anc j (64) _ i ns tead of that, applying the formula lfO|) . 
we can introduce more irrelevant terms in Eqs. (63) and (64) [see the replacement (|52|l ]. 
For instance, we might replace cos(y / 7T0 — nMx/ao) in Eq. (|64a|l with 



^2 D n cos [(2n + l)(v / 7T0 - 7rMx/a )] , 



(65) 



71=0 



where 7? n is a nonuniversal constant. 
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4-3. Symmetries 

Through several considerations, we obtained a field-theory formula for spin operators, 
Eqs. (51) and (64). Using the formula, transformations 1441) - (|46[) and the effective 
Hamiltonians (|41() and H42I1 . we can consider how the symmetries of the spin-1 AF 
chain are represented in the partially-bosonized effective theory framework. 

(i) As we discussed already, a U(l) rotation around the spin z axis, Sj — ► e %1 Sj~ , 
could be realized by 

r}_^e^n, 0^0- JL. (66) 

(ii) The one-site translation Sj — > Sj+i may correspond to 

r)-(x) -> -r]-{x + a ), £,l. R (x) -£l jR (x + a ), 

9(x) -> 9(x + a ) + 0? f cr 3 (x) -> -<r 3 (x + a ) ,„ 7 , 

4>(x) — > 4>(x + a ) — -JttM ' \ /i 3 (x) fi 3 (x + a ) ' [ ' 

These transformations cooperatively change the staggered factor (— I) 3 in front of 
N{x) into 

(iii) The site-parity transformation Sj — > S-j might be regarded as |51| 

ri„{x)^±ii 1 ^{-x), { H ; ^ K 



(68) 




TiKL / cr 3 (:r) -> cr 3 (-x) f Cl -> T&i 



x) ' \ M 3 (x) -> ^ 3 (-x) ' 1 Cb - ±C 



These symmetry operations might tell us which additional terms are allowed to be 
present in formulas (51) and (64): for example, the staggered component N 3 (x) may 
include c 3 ({(iCz. + C,R)d x 7]^} n e m ^ e + h.c), er 3 ^!?^, etc (n: an integer). These terms 
are properly transformed for all the symmetry operations. 

If we integrate out massive fields in the partition function and replace all massive- 
field parts in spin operators with their expectation values, the effective Hamiltonian 
becomes the Gaussian model 1)420 and the reduced spin operators containing only 
massless boson fields are given by 

J 3 (x) ]=d x (l) + CnCOs(V^n(j)-2TTMx/ao] H , 

N + (x) -> e~ i ^ e {fi s ) Di + D 2 cos (Vtocj) - 2nMx/a ) 

(J+(x)) 2 -> e - 2i ^ e C x + C 2 cos (VI^ - 2irMx/a Q ^ + • • 

(N 3 (x)) 2 -► D x + D 2 cos (V4^(f) - 2TrMx/aoJ + • • • , (69) 

where C n and D n are nonuniversal constants, all the Klein factors are naively 
eliminated, and we used (a 3 (x)) 2 ~ const, (r)-rf_) ~ 0(oq 1 ), and (£ 3 £ 3 /) ~ C( a o 
(see |Appendix T3| . This result indicates that the compactification radius of <f> and that 
of 9 are respectively l/yAin and In this bosonized-theory framework, we could 

interpret that the U(l) rotation Sj~ — * e' n S^ , the one-site translation Sj — > Sj+i, and 



Bosonic representation of spin operators 19 

the site-parity operation Sj — > 5_j , respectively, corresponds to 

e^6--^, (70a) 

6(x) -> 6>(x + a ) + Vtt, ^(as) -> 0(sg + a ) - v^M, (706) 

</>(x) -> -0(-2r), 0(x) -> 6(-x). (70c) 

The symmetries of Eq. (70) strongly restrict the emergence of the relevant vertex 
operators in the effective field theory l|42|l : the U(l) symmetry (|70a|) [the translational 
symmetry (f706|l ] forbids e m ^ e [ e *™v^] to exist in the effective Hamiltonian |57|. 
Consequently the TLL state stably remains. 

One can see that the above transformations (70) are very similar to those of the 
effective bosonization theory for spin- 1/2 AF chains (see | Appendix C| |. 

4- 4- Asymptotic behavior of spin correlation functions 

Applying the formulas (51) and (64), let us investigate equal-time spin correlation 
functions in the field-induced TLL phase. The asymptotic forms of equal-time two- 
point functions of J and N are evaluated as 

/ M\ 2 1 

(J 3 (x)J 3 (0)) » — +-(d x <f>(x)d x <t>(0)) 
\a J 7T 



- A" cos UirM-^j ( e -*/EF*<«) e <V5F*(o)} + . . . 



M\ 2 K 1 



ao / 2n 2 x 2 

■A z cos (l-KM^j ( ^) +•••, (71a) 



(J+(x)J-(0)) ^Acos(nM^\ ( e -iV*(*±9)(*) e ^M*±9)<0)) 

x («i/fl(^£/«(o)) + - 



K + K 



(iV 3 (a;)-/V 3 (0)) m 6" cos ^rilf^-j (a 3 (x)<r 3 (0)) 

Jf/2 



cos ( ttM— ) I — ] #„(£)+..., (71c) 



CLQ J \x 

(N+(x)N-(0)) » 6i(^ 3 (s)/i 3 (0))(e-' i ^( a: ) e 4 ^( ) 



£ 2 COS ( 27rM^J ( e -iV47r^(x) e ivW^(0)^ 

c/\ 1/(2K) / i 

— + cos 2ttM— 

x J V a o 



(71d) 
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where and B^'f 1 are nonuniversal constants, K,. is the modified Bessel function, 

(n) ' v ' 

and x ^S> £ c (= vq/tu). In the calculation of Eq. (71), we used Eqs. (|10|l [see the 
comment |48|]. <|A.8I) . I|B.3|) - (|B.5|) . From this result, longitudinal and transverse spin- 
spin correlators are determined as 

(s;ss) « m 2 - A (^) 2 + ^008(2^-) (^) 2K 

/fln\ ( K + 1 )/ 2 

+ T>{ cos (7r(l + M)j) y—j e~ x ^ + ---, (72a) 

/a „- x (K+K- 1 + l)/2 

{S+S„-)«C,co 8 (»A/i)(^) e-'i- 



X 

1/(2*0 



+ V 2 cos(n(l + 2M)j) (^y* +■■■, (726) 

where x = ja 3> £ c , and Cq and T>^ are nonuniversal constants which are 
related to C ( n z) and D { n z '' ] in Eqs. (51) and (64). The first two terms in (S*S§) 
and the second term in (S^~Sq) can also be derived by the NLSM plus Ginzbrug- 
Landau approach [23 12] > but it is difficult to obtain all the other terms within the 
same approach. One can verify that the critical exponent of the incommensurate 
part around wave number p ~ ±2A;_f (Uf = ttM/oq) in (S^Sq), rj z = 2K, and 
that of the staggered part in (S^Sq), rj = 1/(2K), satisfy the famous relation 
f]f] z = 1 |581 1591 12*31 1291 110) . In addition, it is found that the contribution around 
p ~ 7r/a ± kp in (SjSq) and that around p ~ ±kp in (S^Sq) exhibit an exponential 
decay. Comparing Eq. (72) with the spin correlators of a two-leg spin-1/2 ladder in a 
uniform field in Ref. [21] would be instructive (although our calculations in Eqs. (71) 
and (72) are rougher than those in Ref. |2H])- Besides Eq. (72), using the formulas (51) 
and (64), one can calculate various physical quantities (susceptibilities, dynamical 
structure factors, NMR relaxation rates, etc) in the TLL phase |24| . 

In all the calculations of this subsection, we assumed that the three systems H[<f>], 
7i[£ 3 ], and H[n-] are independent of each other. Although small interactions among 
these systems would actually be present, it is expected that their effects in the low- 
energy, long-distance physics are almost negligible and could be absorbed into some 
parameters such as K, f(o), tn, etc. 

5. Some applications of the field-theory representation of spin operators 

In this section, utilizing the derived formulas (51) and (64), we discuss some topics 
for the low-energy physics in/around the uniform-field-induced TLL phase. 

5.1. String order parameter 

As a quantity characterizing the Haldane phase (H < m), there is the nonlocal string 
order parameter which can detect the so-called "hidden AF long-range order" in the 
phase [HOI EH EI - It is defined by 

O a (i,j) = _ lim / Sfexp \m V S« ) Sf) . (73) 



\i— i|— *oo 



/ ; " j 3 
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In the contimious-field-theory framework, we can predict |411 I4U) that the string 
parameter is approximated as follows: 

O a (i,j) ~ (^+ 1 (x l ) f i a+2 (x l )^ +1 (x J )^ +2 (x J )), (74) 

where Xi = iao and Xj = ja^. Actually, the right-hand side is a finite, non-zero value 
in the Haldane phase where (fi a ) ^ 0. 

From the formula (|A.6|1 . the z component of the string parameter is rewritten as 

O z (i,j) ~ (cos( % A : $(x,))cos(v / ^$(^))). (75) 

In the field-induced TLL phase (H > m), the vertex operator exp(±zy / 7r<E>) is 
represented by using fields (</),9,r)-) [see Eqs. (|636|) and Q64afl ]. Hence, it is predicted 
that in the TLL phase, O z (i,j) behaves as 

O z (i,j) ~ (cos(\^tt cj}(xi) — irMi) cos(v / tt '4>{ x j) ~ nMj)) +••• 

f a \ K/2 

~ cos [nM(i - j)} — +■■■. (76) 

\Xi — Xj J 

Namely, O z (i,j) is shown to exhibit power decay in the TLL phase. 
5.2. S\J (2) -invariant perturbations 

In Sees. 15.215.41 we investigate typical perturbations for the critical TLL state. In 
particular, we focus on whether or not the perturbation terms yield a first excitation 
gap, and the symmetries of them. 

In this subsection, we discuss two terms: the bond alternation J^j'^j+l 

(| S | <C 1) and the next-nearest-neighbor (NNN) exchange J2j faSj ' (l-^l <C J), 
which are invariant under the global SU(2) transformation. Because a global U(l) 
symmetry, a part of the SU(2) one, is usually necessary for the realization of the 
TLL i.e., a c = 1 CFT PEHHSj, and (as we mentioned in Sec.Ol it prohibits all 
vertex operators with the dual field 9(x) from emerging in the effective Hamiltonian, 
we expect, without any calculations, that the TLL phase survives even when these 
perturbations are applied. Let us consider the two terms below in more detail. 

From the continuous- field formula l|7jl. we can expect that the bond alternation 
term is approximated as 

^(-i)v^- • s j+1 = £ js(-iy [i (s±sr +1 + h. c ) + s*s* +1 

j 3 

- JSa J dx|(-l) J '[i (j + {x)J-(x + a ) + h.c) + J z (x) J z (x + a ) 
~ (N+(x)N-(x + oq) + h.c) + N z {x)N z {x + a ) 
+- (-J + (x)N-(x + oq) + N + (x)J-(x + ao) + h.c) 

- J z (x)N z (x + a ) + N z (x)J z (x + <z ) j. (77) 

Let us substitute the formulas (51) and (64) into the above result, although such a 
procedure sometimes causes mistakes (see the comment |55|V As a result, we obtain 

/OQ 
dx [K cos(n(V4^0 - 2irMj) + a' n ) 

71=1 
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+ (-1) J A„ cos(n(\/47r0 - 2%Mj) + a n ) 

+ ■■■, (78) 

where An and a.n are nonuniversal constants. Here, we used (<r 3 (x)) 2 ~ const, 
and the operator product expansion (OPE) ^i 3 x £ 3 R ~ cr 3 + • • • ^ |2 EJ I52| . 
and then integrated out all the massive-field parts in the partition function (see 
the comment jSS])- In order to more restrict the form of Eq. J7SJ), we utilize the 
symmetry argument. For the one-site translation and the site-parity operation, the 
bond alternation term changes its sign. Its bosonized form should also have the same 
property; namely, we require the right-hand side in Eq. I|78|) to change the sign for 
<fi(x) — ► <fi(x + ao) — \/ttM and <fi(x) — » — <j){— x) [see Eq. (70)]. Consequently, we set 
\' n = and a n = it/ 2. The resultant form of Eq. (|7H|l is 

JS Sj ■ $3+1 ~> JS I dx^T{-iy\ n sin(n(V4^> - ^Mj)) 

1 J n=l 

+ • • • . (79) 

This bosonized form indicates that in general, the bond alternation is irrelevant in 
the TLL phase due to the staggered factor (— l) 3 and the phase 2nirMj. However, 
when M — 1/2 (half of the saturation), a relevant interaction sin(V47r0) originates 
from the n = 1 term in Eq. I|79|l because of the cancellation of two factors (— l)- 7 and 
2%Mj. The scaling dimension A s of sin(V47r</ ) ) is K (1 < K < 1.5 25 ). At this case 
of M = 1/2, the low-energy physics may be described by a sine-Gordon theory, and an 
infinitesimal bond alternation induces a finite excitation gap and a finite dimerization 
parameter (S^n • <S2n+i — 62 n +i ■ S l 2n+2). The existence of the gap further means that 
the bond alternation brings an M = 1/2 plateau in the uniform magnetization process. 
The scaling argument near a criticality |63| shows that the bond-alternation-induced 
gap Ag and dimerization parameter respectively behave as 

(S 2n ■ S 2n+1 - S 2n+1 ■ S 2n+2 ) ~ -s%n(8)\5\ K IV- K \ (80) 

for a small \5\. Because of the inequality 1/(2 — K) > 1, the gap gradually grows with 
increasing \S\. These predictions for the (small) bond alternation are consistent with 
previous numerical |(S4| and analytical |65) works. 

Next, let us tern to the NNN coupling (J 2 ) term. Through an argument similar 
to that above, we arrive in the following result: 

y^ y J 2 Sj ■ S j+2 ~> J 2 J dx^^-{d x (f>) 2 - ^-j=d x (j> 

OO 

+ ^A„cos(ra(V^-27rMj)) + •••}, (81) 

n=l 

where X n is a nonuniversal constant. We used the symmetry argument: the right-hand 
side in Eq. Ij81|l is invariant under 4>{x) — * <fi(x+ao) — ^/ttM and 4>{x) — + —<f>(—x). Since 
the bosonized NNN coupling does not contain any relevant operators for arbitrary 
magnetization values, we conclude that the TLL phase remains even when a sufficiently 
small NNN exchange perturbation is introduced. The derivative term d x (f> is absorbed 
into the Gaussian part via cf> — > <f> — const J 2 x x and provides a small correction of the 
magnetization M. While the boson quadratic term (d x <fi) 2 makes the velocity v and 
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the TLL parameter K modify. After easy calculations, we can see that when J 2 > 
(J2 < 0), M and K decrease (increase), but v increases (decreases). 

5.3. Axially symmetric terms 

In this subsection, we consider three kinds of U(l)-symmetric perturbation terms: 
the single-ion anisotropy D z ^2j{Sj) 2 (so-called D term), the XXZ type anisotropy 
JA z Y^j SjSj +1 , and the staggered-field Zeeman term along the spin z axis 
—h z J2j(—l) j Sj (\D z \,\h z \ <C Jand|A 2 | <g; 1). As in the cases of the bond alternation 
and the NNN coupling, there is a high possibility that the TLL phase survives as these 
perturbations are added. 

Following the similar argument to that in the last subsection, we can bosonize 
the three terms as follows: 

3 

oo 

+ ^ d n cos(n(V^(f> - 2nMj)) H j, (82a) 

n=l 

JA Z S Z 3 S Z J+1 -> JA Z J dx[^{d x cpf - ^d x </> 

^„cos(n(v / 4^0-27rA/j)) + ---|, (826) 

/OO 
di{(-l) J ^z„ cos(n(\/4^> • 



OO 



+ ^ 



2nMj)) 



n=l 



(82c) 



where d n , S n and z n are nonuniversal constants. For instance, we required the 
bosonization form l|82cll of the staggered-field term to change its sign for the one- 
site translation <f>(x) — > 4>{x + ciq) — y^nM and to be invariant under the site-parity 
transformation <f>{x) — </>(— x). The results (|82 a|l and l|826(l suggest that the D term 
and the XXZ exchange play almost the same roles in the low-energy, long-distance 
physics of the TLL phase. If D z > (D z < 0), M and K decrease (increase), while 
v increases (decreases). As expected, these three terms generally do not destroy the 
TLL state. However, as in the case of the bond alternation term, when the uniform 
magnetization M becomes close to 1/2, the staggered-field term involves a relevant 
term cos(v / 47r^) with A s = K. Therefore, a staggered-field-induced gap A; lz opens 
and (Sp obtains a staggered component at M — 1/2. From the standard scaling 
argument, the gap and the magnetization are shown to behave as 

A hz ~ \h,\W-«\ 

{S z 3 ) « 5 + (-iyf (j ) sgn{h z )\h z \ K /^- K \ (83) 
where / is a function. 
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5.4- Axial-symmetry-breaking terms 

Here, we discuss three axial-symmetry-breaking terms: the D term with the x 
component of spin D x J2j(Sj) 2 , another kind of the single-ion anisotropy (so-called E 
term) ££ J -[(S';) 2 -(Sj) 2 ] = f T,j[( S t) 2 + ( S 7)% and the staggered-field term along 
the spin x axis — h x J2j(~ tyS* (\D X \, \E\, \h x \ <C J). Since the three terms destroy 
the axial U(l) symmetry (9 — > 9 — j/y/n), vertex operators with 9{x) are allowed to be 
present in the effective Hamiltonian. ft is inferred that such vertex operators cause an 
instability of the TLL state, and then a finite excitation gap occurs. Here, note that 
the D x and E terms are invariant under the 7r rotation Sf — ► —Sj (9^9 — ^/n) PD|. 
whereas the h x term obtains a minus sign via the same rotation. Furthermore, the 
7r/2 rotation leaves the E term change the sign. 

Through some calculations, the three perturbation terms are bosonized as 

/oo 
dx( cos(V^0) d n ] cos(n(v / 4^<^ - 2i:Mj)) 

OO 

+ cos{4^6) Y d n ] cos(n(V^r<f> - 2nMj)) 

n=0 

oo 

+ Y d n ] cos(n(%/4^0 - 2irMj)) + •••}, (84a) 
EJ2^ S ] ) 2 - (^j) 2 ] ^ E I dx { cos(%/4^0) x 

3 J 

OO 

^ e„ cos(n(\/4^> - 2nMj)) H | , (846) 

n=0 

- /,,.^; 1 : .v; -► -h x f dx{ cos(V^9) 



oo 

x ^ x n cos(n(V47T(f> - 2nMj)) H |, (84c) 

n=0 

where rf„ , e n and nonuniversal constants. (Since the bosonized form of the 

anisotropic exchange JA X J2j Sj^j+i ^ s tne same type as Eq. 1)84 of . we do not discuss 
it here.) One should note the following properties: (i) cos(v / 7r#) — * — cos(^/n9), 
cos(-\/47r6') — * cos(V4tt9) and cos(4\/7r0) — » cosiAy/nO) for the 7r rotation, and (ii) 
cos(-\/47r6') — * — cos(V4tt0) for the 7r/2 rotation. 

The most relevant operator in both D x and E terms is always cos(\/47r6') with 
A s = 1/K < 2. Thus, the low-energy properties can be explained by a sine- 
Gordon model, and a gap emerges. Supposing and eo are positive, the potential 
cos(^/4tt9) pins the phase field 9 to ±yfir/2 (0 or y/n) modulo 2^/tt for D X ,E > 
(D X ,E < 0). In such a case of 9 — > ±y/ir/2 (— > or y^7r), it is expected that 
(S*J) - (-lp(cos(0F(9)) = 0) and (S$) ~ (-lp(sin(0F6»)) ^ (= 0). From 

this prediction and the scaling argument [A s of e 1 ^ 9 is 1/(4K)], we conclude that 
for a small D x term, the gap A^ x and the transverse component of the spin moment 
increase as follows: 

An ~ {DJW-V, 
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(S x ) ~ 

(4) ~ (-ly^ 1 /^- 4 ) 



(for D a > 0), 



^ 1 ij '^J , (for D x < 0). (85) 

Of course, for a small E term, the similar results hold: we may replace D x to E 
in Eq. (JSSJ. The staggered moment (i.e., a Neel order) along the spin x or y axis 
shows that the D x (or E) term causes the spontaneous breakdown of the one-site 
translational symmetry. Because of K/(2K — 1) < 1 and l/(8K — 4) < 1, the D x - 
or _E-term-induced gap and the transverse moment (Sj' v ) rapidly increases with the 
growth of \D X \ or \E\. 

The staggered-field h x term contains the relevant term cos(y / 7r#) with A s = 
1/(4K) < 2. Therefore, the /i^-induced gap and the staggered magnetization are 
shown to behave as 

A^ ~ \h x r/^\ (Sf) ~ (-iy S g*(h x )\K\^ SK - l l (86) 

For the spin-1/2 AF Heisenberg chain, which low-energy sector is also described by 
a TLL theory (see Introduction), a staggered field also yields a gap and a staggered 
moment. Oshikawa and Affleck |)B] show that in the spin-1/2 case, A/j x ~ l^zl 2 / 3 and 
(Sf) ~ (-I)- 7 'sgn(/i 2 .)|/i 2 .| 1 / 3 . The result JSgJ thus indicates that the growth of both 
the gap and the staggered moment in the spin-1 case is much sharper than that in the 
spin-1/2 case. 



5.5. Magnon decay induced by axial-symmetry-breaking terms 

In this subsection, we briefly mention the magnon-decay processes, which have already 
in some detail discussed in Ref. [201 - I n order to consider such processes, let us go 
back to the effective Hamiltonian (|36fl . where 77+, fjl and fjl_ respectively denote the 
S z = +1, and —1 magnon creation operators. The marginally irrelevant A term, 
omitted in Eq. (|36|) . just contribute to the magnon decay. 

First, we focus on the U(l)-symmetric AF chain (JIJ without any perturbations. 
Since the one-magnon excitations are present only around p = 7r/a , only the decay 
from a magnon to an odd number of magnons is possible. If H is increased so that 
3e+(0) < e_(0) (H > to/2) [3e+(0) < e (0) (H > 2m/ 3)} are satisfied, a magnon fjl 
[ifo] is energetically permitted to decay into three magnons 1)+?)+?)+ via the A term. 
However, since 77+, 77J and ij_ possess different eigenvalues of S z (namely, three kinds 
of one-magnon states are in different sectors of the Hilbert space), this type of the 
decay is forbidden. Indeed, for a U(l) rotation — ► e 11 ' fj- [770] is transformed as 
rj- — > e %1 f]_ [r/o — » ?7o (invariant)], while the product f] + f] + f] + obey a different rotation 
~ ¥ e3l7 V+V+V+- Three kinds of magnons (fj + , fj_) therefore would be 
well-defined quasiparticles in the U(l)-symmetric system (0J. 

On the other hand, when a U(l)-symmetry-breaking term is introduced, there 
is a possibility that the above decay processes are allowed. In the case with the 
D x or E terms, the it- rotation symmetry (S 1 ^ — > 7 = tt), a part of the 

U(l) rotation, survives. For this rotation, f)- and 77+7^+77+ are odd, whereas 770 is 
even. Furthermore, from Eqs. i|38|) an d (|39fl . w e find that for example a sum of 
four terms u i]^(ki)fj + (k2)fi+(~ki)f/ + (— kz) — (fci,2 — * — fa, 2) + h.c" is invariant under 
both the one-site translation and the site-parity operation. (Because the site-parity 
operation (|39|l is the result of the approximation, the requirement of the invariance 
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under this transformation might be too strong.) As a result, the process ff_ — » 77+77+77+ 
is permissible. It is hence inferred that as sufficiently strong D x or E terms are present 
in the system, ryl magnons become ill-defined particles, and we should eliminate the 
magnon fields from the Hamiltonian l|41|) and the field-theory formulas of the spin, 
Eqs. (51) and (64). In the case with the h x staggered-field term, symmetries of the 
7r rotation and the one-site translation are also broken (the two-site translational 
symmetry remains). Therefore, the magnon decay would be more promoted. When 
H is further increased and the 7y+-magnon condensation occurs (H > m), other types 
of the magnon decay are energetically admitted: for instance, 77I — > 77J77J77+77+77+. 

From the simple discussion above, one sees that as H is sufficiently strong, a 
large axial-symmetry-breaking perturbation tends to make the lifetime of massive 
magnons shorten. However, if such a perturbation is small enough, our effective theory 
framework would still be reliable and have the ability to explain various low-energy 
properties of spin-1 AF chains. 

6. Summary 

In this paper, based on the Majorana fermion theory, we have reconsidered the field 
theory description of the spin-1 AF chain J3J|, and derived an explicit field-theory form 
of spin operators in the uniform-field-driven TLL phase in the chain (0J, i-e., Eqs. (51) 
and (64) [the corresponding effective Hamiltonian is Eqs. (|41|l and 1|42|) ]. From the 
formula, we have completely determined the asymptotic forms of spin correlation 
functions (Sec. 14.41) . Furthermore, applying the formula, we have investigated the 
string order parameter and effects of some perturbation terms (the bond alternation, 
the next-nearest interaction, anisotropy terms) in Sec. We have estimated the 
excitation gaps and some physical quantities (staggered moments and the dimerization 
parameter) generated from the perturbations. From Sec. El to Sec. we have often 
argued how symmetries of the spin-1 AF chain are represented in the effective field 
theory world. 

Our results, especially Eqs. (51) and (64), must be useful in analyzing and 
understanding various spin-1 AF chains with magnons condensed (i.e., with a finite 
magnetization) and the extended models of them (e.g., spin-1 AF ladders, spatially 
anisotropic 2D or 3D spin-1 AF systems). The results of Sec. 03 guarantee this 
expectation. We will apply the contents of this paper to other spin-1 AF systems 
in the near future. 

Determining nonuniversal coefficients in Eqs. (51), (64) and l|69(l . especially those 
in front of terms including only massless bosons, is important for more quantitative 
predictions of (quasi) ID spin-1 systems. A powerful way of the determination is to 
accurately evaluate the long-distance behavior of spin correlation functions by means 
of a numerical method such as DMRG and QMC |10| . 
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Appendix A. Abelian bosonization for fermion systems 

Here, we briefly summarize the Abelian bosonization. As mentioned in Sec. [21 in 
(1+1)D case, a massless Dirac (complex) fermion or two species of massless Majorana 
(real) fermions (critical Ising models) are equivalent to a massless bosonic Gaussian 
theory. The former Hamiltonian is written as 



H[ip] = j dx ic (i^^dxipL - il> R d x %l)Rj 



E^M-M)^!^^]. (A.l) 

"3=1,2 

where ip v and Q are respectively the chiral components of the Dirac fermion and the 
real one, and c is the Fermi velocity. These fermions obey anticommutation relations: 

&(y)} = S v ^6 qtg >5(x - y), {M^VAv)} = ° and {Mx)^l(y)} = 
8 ll y5{x — y). The corresponding Hamiltonian of the Gaussian theory is 



2] 



H[4>] = J dx C - [(d x 9) 2 + (d x c^) 2 ] = / dxc[(8 x ^ L ) 2 + {dM 

(A.2) 

where = 4>l + <t>R is the real scalar field, 6 — 4>l — 4>r is the dual field of 
4>, and 4>l(r) is the left (right) moving part of <f>. The fields <f> and 6 satisfy the 
canonical commutation relation [(f>(x),9(y)] — — ^sgn(x — y). Chiral fields 4>l,r obey 
[<f> L (x),<f> R (y)] = and [0ur(x), <f> L/R (y)} = =F|sgn(a; - y). In condensed-matter 
physics, the Hamiltonian (|A. 1|) usually originates from a microscopic system in solids 
(e.g., a lattice system such as the Hubbard chain and the Heisenberg one) after a 
coarse-graining or a renormalization procedure. 

Among these fermion and boson fields, operator identities hold. The fermion 
annihilation (or creation) operators are bosonized as 

ipL = — f= exp(-?v / 47T^L), tpR = — f= exp(i\Z^K<j) R ) , (A. 3) 
\Z2na V27ra 

where kl.r are Klein factors which satisfy {k v ,k V '} = 2i5„.„<, and are necessary for 
the boson vertex (exponential) operators to reproduce the correct anticommutation 
relation between ipL and ip R . The parameter a is a short-distance cut off, which 
depends on details of the microscopic model considered. (Note that it is possible to 
construct another formula without Klein factors, although it requires a modification 
of commutation relations among bosons. See Refs. |57| |HS].) Following Haldane's 
harmonic- fluid approach |f)9l ll()j , one can obtain an alternative bosonized form of 
tpL and if) R : when a real-space fermion field ip(x) in the considering microscopic system 
is approximated as tp ~ e lkpX ip R + e~ lkpX ipL, one may bosonize ipL R as 




& ih F x^ R _ e i{2n+l)(^4,+k F x) e -iV^s _ (A.4) 

v n— 

The quantity fcf is the Fermi wave number in the microscopic system. The n = 
most relevant terms correspond to Eq. (|A.3(1 . The chiral U(l) currents J L =: tp^ipL '■ 
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and Jr =: fp^tpR '■ are written as 

J L = -^=d x <p L , J R = ^=d x (j) R . (A.5) 
V n V 77 

In addition to the correspondences between the fermion and the boson, it is 

known 01 EHl HH E21 EDI that the Ising order and disorder fields a v and \i v can 

be bosonized as [7J 

(T1CT2 ~ sin(-v/7r0), ~ cos(v / 7T'/>), 

ciA*2 ~ cos(v / 7r^), Mi ^ ~ sin(v / 7r^)- (A. 6) 

Following these Abelian bosonization rules, one can bosonize ID interacting Dirac 
fermion systems as well as the free massless fermion (lA.lfl . If the interaction terms 
are all irrelevant in the sense of the renormalization group, the effective Hamiltonian 
at the low-energy limit is still a Gaussian type with the velocity c and coefficients 
of (d x <fi) 2 and (d x 9) 2 renormalized. Conventionally [3], the resultant Hamiltonian is 
written as 



= / dx c - 



g(d x 9) 2 + -(d x( f>) 2 
9 



(A.7) 



where d is the renormalized velocity, and g is called the TLL parameter (the 
Hamiltonian 1A.20 corresponds to a g = 1 theory). When a system is reduced 
to this type at the low-energy limit, we say that the system belongs to the TLL 
universality [1 EJ E| . The Gaussian theory P |2 H E| y ields 



(9 x ^(x)S x </)(0)) = -^i 



,9 9/2 



(^^(o)) „ Sq q , {^) q2/(29) , (A.8) 

where q and q' are an integer (note the comment |5(J|L This result shows that the 
scaling dimensions of d x <j>, e lq ^^ and e lq ^ e are respectively 1, q 2 g/4 and q 2 /(4g). 

Appendix B. Correlation functions of massive fermion theories 

In this Appendix, we evaluate correlation functions of two massive fermion systems, 

M 3 } = J dx % -v (e L d x e L ~ + m»$& 

= j2 e omi(k)m(k), (b.i) 

k 

H[V-]= J dxrl (-^-d 2 x +m + H + 0(d$))r)- 

= J2t-(k)fiUk)v-(k), (B.2) 

k 

where R is the chiral real fermion field, and other fields ?7_ and 770,- are defined in 
Sec. 
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First, we consider the Majorana fermion system IjB.lfl . At T — 0, {%(k)fja(k')) = 
and (fjo{k)fjo(k')) — 8k,k'- Therefore, the two-point function of is calculated as 

<£i(*)£(0)> = ~ £ a^jfe) [ (£o(fc) - + (£o(fc) + fcw «) e ~ 4fc 1 

1 , 1 f A dk kv _.. 

2m 

w -i tfifr/k), (atz»£ c ), (B.3) 

7Tl> 

where we used Eqs. l|3*Tjl. <|3"l|) and (f53f> . Here, A is the ultraviolet cut off, £ c = vo/m, 
and K v is the modified Bessel function [K u [z) ~ \J^~ e ~ z at z > 1). Similarly, one 
can obtain 

{e R {x)e R m = {ti^im, 

(&(x)&(0)) « - i^-K (x/^ c ), (at x » (B.4) 

Z7TW0 

In another system (IB.2|) . the similar relations (rj_(k)rj_(k')) = and 
(fj-(k)ij_(k')) = Sk.k' hold at T = 0. One hence easily finds 

(r ? l( a; )7 7 _(0)) = 0, 

(r ? _(x)^(0))=5(x)~^, (B.5) 
a 

where a is the short-distance cut off. 

In addition to these results, one can of course compute any correlation functions 
of the systems (|B.1|I and (|B.2|) . using Wick's theorem, etc. 

Appendix C. Symmetries of the spin-1/2 XXZ chain 

The spin-1/2 XXZ chain has the same global symmetries as those of the spin-1 
AF chain Q): the U(l) rotation around the spin z axis, the one-site translation, and 
the site-parity transformation. In the Abelian bosonization framework (see Eqs. (J2J 
and 121), these three symmetries could be realized by the following transformations of 
boson fields and 6 [Z1E2- 

(i) The U(l) rotation Sj~ — > e 11 S~!f corresponds to 

9^9 + ^=. 

(ii) The one-site translation Sj — > Sj+i corresponds to 

r 



(C.l) 



0(x) -> 0(a; + o ) + I M . , ) 

0(ar) -► 0(z + ao) + v^ 7 . (C.2) 
(iii) The site-parity transformation Sj — > S—j corresponds to 

0(a:) -» (C.3) 
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